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This paper presents the analytical and numerical analysis of fluctuational dynamics of ac hysteretic 
SQUID. It has been demonstrated, that the most important parameter for the improvement of the 
hysteretic SQUID sensitivity is the ratio between the pump and the characteristic frequencies but 
not their absolute values. The resonant behavior of signal-to-noise ratio as function of the pump 
frequency has been observed and the narrow area for the optimal pump frequency range is indicated. 



It is well known |E]],|§|,0] that the slope of 
a signal characteristic of hysteretic ac SQUID 
h = dV ac /d<P increases with increasing the working 
(pump) frequency and therefore output energy sen- 
sitivity becomes considerably better. Noise prop- 
erties of ac SQUIDs were intensively investigated 
both for cases f2 <C 1 and £1 ^> 1 (where SI = uj/uj c , 
lj is the pump frequency and lo c is the characteris- 
tic frequency of Josephson junction (JJ)) H,||,||. 
It has been demonstrated that increasing of fre- 
quency without violation of working capacity of 
SQUID is possible nearly up to the characteristic 
frequency of J J M , which may be of the order of 1 
THz for J J based on low T c materials and 10 THz - 
for high T c ones. Besides, it was reported ^ about 
ac SQUID with f2 = 0.1, demonstrated high sen- 
sitivity about 4 • 10~ 32 J/Hz. So, elaboration and 
fabrication of a single junction SQUID working in 
microwave frequency range looks very promising, 
and usual opinion is that increasing of the pump 
frequency may allow to reach, in principle, such 
high characteristics as dc SQUIDs have. 

However, at the present time, there is still no 
at least qualitative understanding at which pump 
frequency a microwave SQUID should operate to 
reach the maximal sensitivity, since due to mathe- 
matical problems no investigations were performed 
in the most interesting pump frequency range u ~ 
lj c . As follows from recent results by Cheska ||, 
there should not be such optimal pump frequency 
for a nonhysteretic SQUID, which is not surprising 
since that is a (quasi) linear system. The ac hys- 
teretic SQUID, is, however, strongly nonlinear sys- 
tem, and, recently, in qualitatively similar model 
system, the resonant behavior of signal-to- noise ra- 
tio as function of frequency of driving signal has 
been observed |To| . 

The present paper is devoted to studying the 
nonlinear fluctuational dynamics of a hysteretic ac 
SQUID in the pump frequency range around oj c . 
The effect of suppression of noise by strong peri- 
odic signal JT3] has for the first time been observed 
in the model, describing the real electronic device: 
it has been demonstrated that there exists a cer- 
tain optimal pump frequency range at which the 
hysteretic SQUID will operate with minimal noise- 



induced error and that there is an optimal driving 
amplitude at which the voltage-flux characteristic 
is the most close to the noiseless case. 

Theory of a single junction SQUID is already 
quite well developed (see, e.g., PJ-H and references 
therein). Necessary requirements for SQUIDs op- 
erating in hysteretic mode are the following: i > 1 
to provide hysteresis and maximize nonlincarity 
(where £ = L/Lq, L is the inductance of the ring, 
Lq = •Pq/(2ttI c ), I c is the critical current of the JJ, 
<Pq is the flux quantum) and (3 = 2ttI c R 2 n C /$o <C 
1 (the SQUID should be overdamped) to prevent 
from stochasticity when switching between states 
with different trapped fluxes in the SQUID ring 
(here R^ 1 = Gn is the normal conductivity of the 
J J, C is the capacitance). 

The fluctuational dynamics of a flux in a SQUID 
ring coupled to a resonator may be described by 
the following equations: 

<p + £sw.(p + ip - ip m - a s ip(t) = <p F (t), (1) 
— *V H — ~4' + i> = ctJip - ip m ) + asm(uj t). (2) 
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Here <p = ——— , <P is the trapped flux, ip rn — — — — , 
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<P m is the measured flux and ip(t) is the pump 
signal from the resonator, time is normalized to 
the characteristic frequency of the SQUID lo s = 
lo c /£ — R N /L, ijj c — 2ttR n I c /^ - characteris- 
tic frequency of J J, r) r is the damping of the res- 
onator, a s and a r are coupling coefficients of the 
SQUID ring and the resonator, respectively, and 
loq = Lj r /oj s is the dimensionless pump frequency. 

Let us take into account only internal thermal 
fluctuations in the SQUID ring. In this case the 
noise source <fF(t) may be represented by the white 
Gaussian noise: 

(<p F (t))=0, (<p F (t)<p F (t + T)) =2 7 «(r), (3) 

where 7 = 2irksT /(<!>() I c ) is the dimensionless 
noise intensity, T is the temperature and ks is the 
Boltzmann constant. 

Let us first consider analytically and numeri- 
cally the fluctuational flux dynamics of a SQUID 



ring, supposing that the signal from the resonator 
is known: a s ip(t) = Asin(u>ot + £), where A is the 
driving amplitude and £ is the initial phase. In this 
case the required statistical characteristics may be 
computed on the basis of the transitional proba- 
bility density W(cp, t). It is well known, that the 
Fokker-Planck equation (FPE) for the probability 
density W(ip, t) corresponds to Eq. (pi) for the flux: 
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u(ip) = -I cos(^) + (ip- ipe) 2 /2 (5) 

is the dimensionless potential profile, ip e = ip rn + 
Asm(uj()t + £) is the external flux. The initial and 
boundary conditions for Eq. (Q) are as follows: 

W{ip, 0) = 5{ip - ip ) and G(±oo, t) = 0. 

The quantity, that reflects the noise induced 
transition process from the metastable state is a 

probability P(t) = / W(ip, t)dip that transition 

will not occur at the moment of time t, that we 
will call the survival probability, where d is some 
boundary point, usually a potential barrier top. 

To describe analytically the survival probabil- 
ity one can recourse to the adiabatic approxima- 
tion, that for the case of a weak periodic driving 
has been used in the context of stochastic reso- 
nance ]li"| ] , |l^ ] . It has been demonstrated in [H| 
(see also ]14|) that for the case of a metastable 
potential and strong periodic driving the evolu- 
tion of survival probability in the frequency range 
< luq < 0.5 may be well described by a modified 
adiabatic approximation, which allows to extend 
the usual analysis to arbitrary driving amplitudes 
and noise intensities. For simplicity, let us consider 
the case for ip m = tt (bistable potential), then the 
survival probability may be written in the form: 

P(<Po,0) ' '"' 



P(tp Qi t) = g(t) 



j(t) = cxp ^ - 



o g(t')T q (ip ,t') 
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where T p (tpo, t') is the exact mean decay time of the 
left metastable state obtained for the correspond- 
ing time-constant potential 15 

t> ) = — t | / e"M/7< / e- u ^dxdy + (7) 



u(y)/~fl dy / e -u{x)hl^ 



Here c > d is the coordinate of the right potential 
minimum at the moment when sin(a>o£ + = 1 
and the left minimum disappears. The time scale 
T q is the decay time of the right metastable state 
and is calculated from the same formula (Q) where 
in the potential (0) the phase of the driving signal 
is changed from to tt. Note that with respect 
to the usual adiabatic analysis the approximate 
Kramers' time has been substituted by the exact 
one (]?]) and a surprisingly good agreement of this 
approximate expression with the computer simula- 
tion results has been found in rather broad range 
of parameters, that is seen in Fig. 1. 




Fig. 1. The evolution of survival probability; dots - 
computer simulations, solid line - modified adiabatic 
approximation. Inset: the mean transition time; dots 
- computer simulations, solid line - modified adiabatic 
approximation; 7 = 0.04; 0.1; 0.2 from top to bottom. 

In Fig. 1 the probability evolution is presented 
for 7 = 0.1, A — 3, £ = and ip m — n. The adi- 
abatic approximation is drawn by solid line, while 
the results of numerical solution of Eq. (|4|) by dots. 
One can see, that the coincidence is good enough. 
Even better agreement between the modified adi- 
abatic approximation (^), (^) and the computer 
simulation results may be observed for the mean 
decay time of the left metastable state (inset of Fig. 
1). Since the working conditions for the considered 
potential are such that during half of the driving 
period T w the survival probability monotonically 
changes from unity to almost zero, we can define 
the mean decay time in time-periodic potential as: 



Tp(Wo) 



J^ /2 [P(t)-P(0)]dt 
P(T u /2) - P(0) 



(8) 



The definition (||) is analogous to widely used defi- 
nition of integral relaxation time for diffusion in 
time-constant potentials (see [O and references 
therein). As it is seen from the inset of Fig. 1, 
in the range of validity of adiabatic approximation 
the mean decay time (0) computed using formu- 
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las (||),(0) gives perfect coincidence with the re- 
sults of computer simulation. It is intriguing to 
see that there is a range of driving frequencies 
0.1 < ujq < 0.3, where the mean decay time is al- 
most insensitive to the noise intensity: the curves 
for different noise intensities actually coincide that 
constitutes the effect of suppression of noise by 
a strong periodic signal. This effect may also be 
demonstrated by the analysis of the signal-to-noise 
ratio R. In accordance with [021 we denote R as: 



R = 



1 



lim 



S(u)du, (9) 



ojo — Aw 



where S(u) = f_°° e~ %UT K[t]cIt is the spectral 
density, Sn(^o) is noisy pedestal at the driving 
frequency lo and K[t] is the correlation function: 
K[t] = {{(f{t + T)ip(t))), where the inner brackets 
denote the ensemble average and outer brackets in- 
dicate the average over initial phase £. In order to 
obtain K[t] the Eq. ([!]) has been solved numeri- 
cally for A = 3, I = 3, ip m = 7r. 
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Fig. 2. The signal-to-noise ratio as function of driving 
frequency. 

Let us plot the signal-to-noise ratio (SNR) as 
function of driving frequency luq. From Fig. 2 
one can see, that R as function of ojq has strongly 
pronounced maximum. The location of this max- 
imum lies below the cut-off frequency (loq = 1 or 
lu s = Rn/L in dimensional notations) and depends 
on the driving amplitude A and noise intensity. 
The existence of optimal driving frequency may 
be explained in the following way: with increase 
of frequency noise-induced escapes occur at lower 
barrier, and the maximum of SNR corresponds to 
the lowest barrier, that is closest to the dynamical 
case, where the transition occur when the barrier 
completely disappears. With further increase of 
frequency, there is not enough time for transition 
in the absence of noise and SNR drops, while the 
mean transition time rises E0| , |l3]] . 



One can solve numerically the system of equa- 
tions (0) and (||) and compute ac current-voltage 
characteristic and voltage-flux characteristic. 
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Fig. 3. The ac current-voltage characteristic for 
ujo = 0.01; solid line -7 = 0, dashed line - 7 = 0.01, 
diamonds - 7 = 0.03. 
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Fig. 4. The voltage-flux characteristic for u)q = 0.01; 
solid line -7 = 0, dashed line - 7 = 0.01, diamonds - 
7 = 0.03. 

In Fig. 3 the ac current-voltage characteristic 



0;7r/2;7r) is presented 



for ljq = 0.01 that corresponds to relatively low 
frequency case and for the following parameters: 
I = 3, r\ r = 0.01, a s a r /r] r = 2. As it is seen, for 
the case 7 = the plateaus are nearly horizontal. 
For 7 = 0.01 and 7 = 0.03 the plateaus have a tilt 
and lie below the curves for 7 = 0. The voltage- 
flux characteristics for the same parameters and 
a = 0.47 is presented in Fig. 4. It is seen that 
with increase of noise intensity 7 the correspond- 
ing curves move down and their edges are rounded. 
Low frequency interference, as well as noise of the 
resonator and amplifiers (that were not taken into 
account in the present paper) will make the situa- 
tion only worse and the quantity of splitting of the 
curves for different 7 will give the error of the mea- 
sured flux ip rn . Fortunately, with the increase of 
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the pump frequency the picture changes radically: 
already at ujq = 0.1 the plateaus for 7 = 0.01 cross 
the plateaus for 7 = at certain points. With 
approaching the maximum of signal-to-noise ratio 
uiq k 0.3 (Fig. 2) the crossing points move to the 
middle of the plateau and for different 7 become 
more close to each other, see Fig. 5. 
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Fig. 5. The ac current-voltage characteristic for 
ujo = 0.3; solid line -7 = 0, dashed line - 7 = 0.01, 
diamonds - 7 = 0.03. 




Fig. 6. The voltage-flux characteristic for u>o = 0.3; 
solid line -7 = 0, dashed line - 7 = 0.01, diamonds - 
7 = 0.03. 

It is intuitively obvious that to reach the max- 
imal sensitivity of the SQUID, the voltage-flux 
characteristic should be computed for the pump 
amplitude a, corresponding to the crossing point 
of plateaus for (p m = tt/2 (see analogous consider- 
ation for dc SQUIDs in [jl6|). In Fig. 6 the voltage- 
flux characteristic is plotted for ujq = 0.3, a = 0.51 
and 7 = 0; 0.01; 0.03. It is seen, that there is a 
broad range of cp m , 0.7 < ip m < 2.2, where curves 
for different 7 coincide and, therefore, noise actu- 
ally do not affect the SQUID in this range of pa- 
rameters. If, however, the measured flux lie outside 
the range 0.7 < ip m < 2.2, some known quantity of 
magnetic flux ip n may be added to shift the work- 
ing point to this region, precise measurements may 



be performed, and ip n may be extracted to get in- 
formation about the original measured flux. 

Up to date tendency in ac SQUID develop- 
ment is increasing the pump frequency up to mi- 
crowave range (~ 10GHz) in order to reduce the 
noise || with respect to measured voltage and 
thereby to increase the sensitivity of devices of 
this kind. However, as follows from the above pre- 
sented analysis, the main optimizational parame- 
ter of the microwave hysteretic SQUID is the ratio 
between the pump and the characteristic frequen- 
cies ujq — u> r £/u) c , but not their absolute values; 
it must be in the range u>o « 0.2 — 0.4 in order to 
get maximal sensitivity of the SQUID. Also during 
measurements, very important parameter is the 
amplitude of ac driving. It can be chosen from 
ac current-voltage characteristic of the SQUID as 
the point, where the plateau for <P m = <l>o/4 for 
a given temperature fc^T ^ crosses the plateau 
for ksT = 0; the calibration can be done once by 
cooling the SQUID to a low temperature. 
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